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Abstract 

We report the first study of the dynamics of point defects, mono and di- 
vacancies, in a confined 2-D cohoidal crystal in real space and time using 
digital video microscopy. The defects are introduced by manipulating indi- 
vidual particles with optical tweezers. The diffusion rates are measured to be 
Dmono/cL^ — 3.27±0.03Hz for mono- vacancies and Ddi/a^ = 3.71±0.03Hz for 
di-vacancies. The elementary diffusion processes are identified and it is found 
that the diffusion of di-vacancies is enhanced by a dislocation dissociation- 
recombination mechanism. Furthermore, the defects do not follow a simple 
random walk but their hopping exhibits memory effects, due to the reduced 
symmetry (compared to the triangular lattice) of their stable configurations, 
and the slow relaxation rates of the lattice modes. 
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Uniform colloidal spheres suspended in a solvent, under appropriate conditions can self- 
assemble into ordered crystalline structures Using these crystals as a model system, 
a great variety of problems in materials science, physical chemistry and condensed matter 
physics have been investigated during the past two decades. However, the properties of 
simple structural defects, such as vacancies and interstitials, dislocations, etc., ubiquitous in 
real solids, have received little, if any, attention. Considering the unique possibility of real 
space and real time imaging, colloidal systems provide a great potential for research in this 
direction. 

Recently, it has been demonstrated that point defects can be created in two- 
dimensional colloidal crystals by manipulating individual particles with optical tweezers. 
Direct imaging of these defects verified that their stable configurations have lower symmetry 
than the underlying triangular lattice, as predicted by numerical simulations for a number 
of two-dimensional systems . The same study revealed that point defects can dissociate 
into pairs of well-separated dislocations, a topological excitation permitted only in two but 
not three dimensions. In this Letter we report the first study of the dynamics of mono and 
di-vacancies in two-dimensional colloidal crystals. We discover the first evidence that topo- 
logical excitations into dislocation pairs enhance the diffusion of di-vacancies. Moreover, the 
hopping of the defects does not follow a pure random walk, but exhibits surprising memory 
effects. 

The colloidal crystals are prepared by 0.360 fim diameter negatively charged polystyrene- 
sulfate micro spheres (Duke Scientific, #5036), suspended in highly deionized water. The 
particles crystallize due to strong electrostatic interaction. A single layer of particles is 
confined between two fused silica substrates separated by ~2 /im, resulting in a single layer 
colloidal crystal with a lattice constant a = l.l/xm. The experimental details can be found 
in P . Point defects are created by trapping a particle with optical tweezers and dragging it 
away from its lattice site. The dynamics of such single, isolated point defects are recorded 
on videotape, using a CCD video-camera. Each individual frame (l/60th of a second) 
was acquired on a PC and processed using a particle-tracking algorithm 0. A triangulation 
algorithm is used to identify the mis-coordinated particles and characterize the configuration 
of the defect. 

We tracked a number of mono and di-vacancies for long enough times (up to 40 seconds) 
to obtain reliable measurements of their diffusion constants. Typical trajectories are shown 
in Fig.^ The position of the defect x{t) is defined as the average of the positions of the 5-fold 
coordinated particles around the core of the defect. x{t) is determined with a resolution 
of ~ 0.1a. Thus one can identify the defect in a particular lattice site or in-between two 
neighboring sites. To estimate the diffusion constants, we take the average of the squared 
displacements < Ax'^{6t) > as a function of the time separation 6t, running through all the 
available trajectories. A linear relationship < Ax"^ > vs. 6t is obtained for the range 1-5 
seconds (Fig.|^A). The data in this range are weighted by their errors and fitted to a zero 
y-intercept linear relation, whose slope gives the respective diffusion constants. The results 
for mono and di-vacancies are: Dmano/a^ = 3.27 ± 0.03Hz, Ddi/a^ = 3.71 ± 0.03Hz 0. 

The diffusion constants measured above give the combined contribution of all the possible 
diffusion pathways. In this Letter, we reveal the role of topological excitations on the 
diffusion of point defects. As was shown in 0, the possible stable configurations of mono 
and di-vacancies are not just the "well-identified" point defect configurations but there 
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are also configurations where the point defects appear as separated dislocation pairs. The 
defects are observed to dissociate into two dislocations that glide from site to site and finally 
recombine at a new position that can be several lattice constants away from the original one. 
The individual dislocations move faster than the hopping of mono or di-vacancies and this 
topological excitation has a lifetime Tq of a fraction of a second. Such a process is a distinct 
diffusion pathway, giving a contribution to the diffusion constant Dpair ~ TcT^^^? where P is 
the mean square separation of the dislocation pair. 

Mono- vacancies spend considerably less time as a dislocation pair than di-vacancies and 
the aforementioned diffusion mechanism is not important in their case. One can clearly 
recognize the different behavior exhibited by mono and di-vacancies, as illustrated in Fig.|^, 
where the separation of the dislocations comprising the defect is shown as a function of time. 
The mono-vacancy is observed primarily to fluctuate between configurations in which the 
dislocations are bound, with very rare excursions into configurations where the dislocations 
dissociate. In contrast, the di- vacancy frequently spends long intervals (a few seconds) as a 
dislocation pair, with the two dislocations dissociating and recombining several times during 
this period 

The combined diffusion constant for di-vacancies is Ddi = pDpair + (1 —p)D%, where p ~ 
0.25 is the fraction of time excited into a dislocation pair, and Dpair-,D'^^ are the contributions 
from the dissociation-recombination mechanism and from the "bare" di- vacancy hopping. In 
order to estimate the diffusion constants for the two different mechanisms, we link into two 
separate trajectories the intervals spent as a dislocation pair and as a "bare" di- vacancy. To 
avoid systematically picking up hoping events, the shortest interval was 1 second (larger than 
the time between jumps). Carrying out the same procedure as for the original trajectories, 
the two diffusion constants are found to be Dpair/ cP' = 5.90 ± O.OOHz and D^Ja? = 2.97 ± 
0.03Hz (Fig.QB). These numbers combine to give the measured value Dcn/a? = 3.70Hz. 

We now turn our attention to the details of the underlying stochastic hopping processes. 
Naively, one would expect the trajectories of the defects to be pure random walks on a 
triangular lattice, i.e. when hopping, the system chooses with equal probability between 
the six neighbor sites and consecutive events are uncorrelated. However, a careful analysis 
reveals a peculiar feature: their motion looks nothing like a 2D random walk, since typically 
it is restrained along particular directions . The existence of easy directions for diffusion 
is related to the low symmetry of the stable configurations of the relaxed lattice around the 
defect PJ^. Therefore, the probability for hopping towards the six nearest neighbors is not 
the same for all of them, but is higher along particular lines of symmetry, Fig.^. For di- 
vacancies excited into a pair of separated dislocations, the easy direction coincides with the 
gliding direction of the two dislocations. The easy direction for diffusion remains unchanged 
for up to a few seconds. 

The existence of easy directions for hopping makes the trajectories of the defects quasi- 
one-dimensional, at least for a few consecutive steps. This effect is demonstrated by measur- 
ing the probability Po{n) for returning to the original position after the defect has executed 
n steps. It is found that the measured Po{n) for small n resembles the one for random walk 
on a 1-D chain, rather that for random walk on a 2-D triangular lattice. This is illustrated 
in Fig.|^: Po{n) clearly deviates from the expected 2-D behavior, following the characteristic 
oscillations of the 1-D random walk of period=2 steps. 

A second surprising peculiarity involves memory during site-to-site hopping. The defects 
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appear to spend some time around a lattice site, attempting to hop towards neighboring 
sites by changing into a spht intermediate configuration between the two lattice points. 
Depending on the symmetry of the relaxed lattice in the original configuration (2 or 3-fold), 
the defect is seen to attempt hopping towards either two or three of the six neighbors. Oddly, 
the defect typically hops back to the original site, finally getting to the next site only after 
several attempts. The back-hopping to forward-hopping ratio is found to be ~ 3 : 1. 

This memory effect originates from the slow relaxation rates of distotrions in the over- 
damped colloidal crystal: the core region of the defect can rearrange rapidly (~ lOifz), 
however the rest of the lattice around the defect cannot follow at the same rate (e.g. distor- 
tions with wavelength 10a relax at ~ IHz). As a vacancy changes from a configuration at 
one site to a split intermediate configuration between two sites, the rest of the lattice does 
not immediately relax, and seems to "push" the vacancy back to its initial site [|1^ . 



In summary, the dynamics of point defects in colloidal crystals are studied in real space 
and time. Di- vacancies in particular are observed to dissociate into two dislocations that can 
separate several lattice spacings before recombination. This distinct mechanism enhances 
the diffusion rate of di- vacancies compared to their "bare" one. The reduced symmetry of 
the defects' stable configurations results in easy directions for diffusion. We have also found 
that the site-to-site hopping seems to exhibit memory effects accounted for by the slow 
relaxation of the short wavelength modes of the crystal. We acknowledge helpful discussions 
with Prof. S.C Ying and Prof. D.A. Weitz. This work was supported by NSF Grant No. 
DMR-9804083, the Petroleum Research Fund, the Research Corporation, and the A. P. Sloan 
Foundation. 
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FIGURES 



FIG. 1. Trajectories of defects: (Top) mono, (Bottom) di-vacancy. Neighboring clumps of dots 
are separated by roughly half a lattice spacing and correspond to the defect being in a configuration 
near a lattice point or in a split configuration between two lattice points. A snapshot of the lattice 
is super-imposed a few seconds after the defect diffused away from the region shown. The positions 
of the defect are displaced with respect to the lattice seconds later due to slow, long wavelength, 
thermally generated distortions. 

FIG. 2. Top: Mean square displacements of mono and di-vacancies. The average of Ax^ as 
function of At was measured from about 2-minute long mono and di-vacancy trajectories. The 
straight lines are linear fits for At 1-5 seconds. Bottom: Ax^ separately for "bare" di-vacancies 
and di-vacancies excited into dislocation pairs. Fits are for At 1-3 seconds. 

FIG. 3. Separation r{t) of the dislocations comprising a mono and a di-vacancy (top and bottom 
respectively). r{t) is measured as the standard deviation of positions of the 5- fold coordinated 
particles in the defect core. The positions of the minima in the dislocation pair potential |Q] are 
indicated by the dotted horizontal lines. The first two correspond to the split and the crushed 
configurations, while above the third one (emphasized with a dashed line) the dislocation pair is 
dissociated. Note that in the mono-vacancy case a large fraction of the points shown (between SV 
and V2) corresponds to the V3 configuration, which is a bound dislocation triplet. In the case of 
di-vacancy, the arrows indicate intervals in which the defect appears as a dislocation pair. The two 
dislocations repeatedly dissociate and recombine during this period. 

FIG. 4. Configurations of mono-vacancy, from Top: three-fold symmetric V3, bottom: 
crushed V2a- The 7-fold coordinated particles are more probable to hop towards the vacant site, 
resulting in particular easy directions for diffusion (indicated with arrows). 

FIG. 5. Po{n) measured from 40 seconds of data or roughly 60 hopping events. Clearly, the 
observed behavior for both type of defects resembles more the one for 1-D random walk, rather 
the one for 2-D. Po{n) exhibits oscillations of period 2 steps and the average trend lies between the 
1-D and 2-D curves. 
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